BABAR-CONF-06/013 
SLAC-PUB-11978 
July 2006 



A measurement of CP-violation parameters in B°B° mixing using 
partially reconstructed D*~£ + f£ events at BABAR 



The BABAR Collaboration 
February 7, 2008 



Abstract 

CP violation in B°B° mixing is characterized by the value of the parameter \q/p\ being different 
from 1, and the Standard Model predicts this difference to be smaller than 10 -3 . We present 
a measurement of this parameter using a partial reconstruction of one of the B mesons in the 
semileptonic channel D*~£ + V£, where only the hard lepton and the soft pion from the D*~ — > D°ir~ 
decay are reconstructed. The flavor of the other B is determined by means of lepton tagging. The 
determination of \q/p\ is then performed with a fit to the proper time difference of the two B decays. 
We use a luminosity of 200.8 fb™ 1 , collected at the T(4S) resonance by the BABAR detector at the 
PEP-II asymmetrical-energy e + e~ collider, in the period 1999-2004. We obtain the preliminary 
result: 

\q/p\ - 1 = (6.5 ± 3.4(stat.) ± 2.0(syst.)) • 10~ 3 
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1 INTRODUCTION 



Although CP violation has been established both in the interference between decay and mixing and 
in direct B decay, CP violation in the mixing alone has up to now eluded experimental observation. 

In the standard mixing formalism, the effective Hamiltonian is expressed as the sum of a mass 
and a decay matrix (H = M — i/2 T) and the B mass eigenstates are connected to the flavor 
eigenstates by: 

\B L }=p\B°}+q\B } 

\B H )=p\B°)-q\B°). (1) 
The absolute value of the ratio q/p can be written in this notation as: 





2 




1 






p 







'M* 2 -i/2 r* 12 



m 12 - i/2 r 



12 



2 



so it is exactly equal to 1 in a CP-conserving scenario (where the off-diagonal elements of the mass 
and decay matrices are real, that is M\i = M* 2 , T12 = T^ 2 ), while it differs by a small quantity 
if CP is violated in mixing. The current theoretical predictions on this quantity in the Standard 
Model (SM) are 2 • 1CT 4 < \q/p\ - 1 < 6 • 10" 4 [BEj. 

However, recent theoretical publications [2| have pointed out that in some New Physics (NP) 
general scenarios, the predictions for this quantity can be significantly different with respect to the 
SM. Making only the assumptions that the CKM 4 is a 3 x 3 unitary matrix, and the tree-level 
processes are dominated by the SM, the matrix element M\i in this scenario can be related to the 
SM one by the formula: 

Mf 2 p = r 2 d e 2ie *Mi™ (3) 

where r& and Od are general New Physics amplitude and phase, while Ti2 is not modified. The CP 
asymmetry in mixing can be written as: 

T(B° -> i + X) - T{B° M) _ 
SL ~ T(B° Z i+X) + T(B° Z £-X) ~ 

* 2(1 -\ q / P \) = -Re ^—j _ r+ lm^— j — (4) 

If the New Physics phase is significantly different from 0, the real part of Ti2/Mi2 can become 
dominant, enhancing the asymmetry up to an order of magnitude. Together with the measurements 
of CKM parameters, an accurate determination of \q/p\ is therefore an additional constraint on 
New Physics models. 

Measurements of \q/p\ at the i?-factories are performed both using inclusive dilepton events 
E] or B mesons fully reconstructed into flavor or CP eigenstates 0. The most precise result 
obtained up to now from the inclusive dilepton method is \q/p\ — 1 = (—0.8 ± 2.7 ± 1.9) • 10~ 3 [5], 
where the first uncertainty is statistical, the second systematic. Analogous techniques have also 
been used by the D0 experiment at the Tevatron |Hj: the measured parameter in this analysis is not 
a determination of \q/p\ alone, but also involves contributions from the CP violation in B s mixing 
[SI, and this has allowed to set constraints also on New Physics in the B s sector. 

In this analysis we exploit the partial reconstruction of B° — ► D*~l + u^ events. Though the 
total statistics is not as high as in the dilepton case, we can keep the charged B background at 
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a lower level, while selecting a greater number of events than in an analysis requiring the full 
reconstruction of a hadronic or semileptonic decay as a tag for the B flavor. At the same time, 
since the reconstructed and tag side are well defined, a procedure to determine particle detection 
asymmetries from data (explained in Sec. QJ) can be carried out. In this way, we do not need 
dedicated control samples to determine the asymmetry induced by the experimental cuts. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetrical- 
energy e + e~ storage ring in the period 1999-2004; they correspond to an integrated luminosity of 
200.8 fb _1 (i.e. about 110 million B°B° pairs) taken at the mass of the T(4S) resonance (y/s = 
10.58 GeV), plus 21.6 fb" 1 taken about 40 MeV below. Samples of simulated T(45) -> B°B° and 
B + B~ events are used to estimate efficiencies, study background and detector asymmetries. 

The BABAR detector is described in detail elsewhere ^Uj. Tracking of charged particles is pro- 
vided by a five-layer silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). Vertices and 
soft pion tracks are mainly reconstructed using information from the SVT. Cherenkov radiation 
detected in a ring-imaging detector (DIRC) is used for particle identification. An electromag- 
netic calorimeter (EMC), which consists of 6580 thallium-doped Csl crystals, is used to measure 
electron energies. These systems are mounted inside a 1.5 T solenoidal superconducting magnet. 
The flux return of the magnet (IFR) is equipped with Resistive Plate Chambers providing muon 
identification. 

3 SELECTION METHOD AND SAMPLE COMPOSITION 

The decay rates for neutral B mesons can be calculated theoretically as a function of At, the 
time difference between the decays of the two B° mesons, taking into account both the time 
evolution of the mass eigenstates and the fact that they are produced coherently from the T(4S) 
resonance. We approximate At ~ Az/(((3^f)c), where Az is the measured distance between the 
decay vertices projected along the beam direction [z axis) and (f3j) is the average boost of the 
T(4:S) in the laboratory frame. The method used to determine \q/p\ is a two-dimensional fit to the 
set of variables (At, a^t) with a binned extended maximum likelihood approach, where c^t is the 
per-event error calculated on At. 

To discriminate events with a BB pair from events which originate from the production of light 
quarks, we require the ratio of the second to the zeroth Fox- Wolfram moment R2 < 0.5. The 
total number of charged tracks in the event is required to be greater than 4, to discard lepton 
pair production. We also calculate the invariant mass of the two highest-momentum leptons in 
the event and we apply a veto on the regions Mu < 0.35 GeV/c 2 (converted photon rejection) and 
3.07 < M u < 3.14 GeV/c 2 (J/ip rejection). 

Different techniques are then used to identify the two B mesons in an event. For the first 
(B rec ) we select B° — > D*~£ + V£ events with partial reconstruction of the decay D*~ —* D ir~, 
using only the charged lepton from the B° decay and the soft pion from the D*~ decay. 

The D° decay is not reconstructed, resulting in high selection efficiency. Electrons are identified 
using dE/dx measurements in the DCH, the DIRC information and the ratio between the energy 
deposited in the EMC crystals and the measured momentum in the DCH. Muons are identified 
using the dE/dx measurements in the DCH, the DIRC information and the number of hits in the 
IFR. Due to the limited phase space available in the D* decay, the tt s is emitted within a cone with 
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a half-opening-angle of approximately one radian and centered about the D* motion direction |12j . 
We approximate the direction of the D* to be that of the ir s and estimate the energy Ed* of the 
D* as a linear function of the energy of the ir s , with parameters taken from the simulation. We 
define the square of the missing neutrino mass as: 



where all quantities are defined in the T(4S) frame. We neglect the momentum of the B° (approx- 
imately 0.34 GeV/c), and identify the B° energy with half the total energy of the events [yfs/2). 
Ep and pp are the energy and momentum vector of the lepton and pD* is the estimated momentum 
vector of the D*. The distribution of is peaked for signal events, while it is spread over a wide 
range for background events. 

We determine the B° decay point from a vertex fit of the t and tt s tracks, constrained to the 
beam-spot position in the plane perpendicular to the beam axis (the x-y plane). The beam spot 
position and size are determined on a run- by-run basis using two-prong events ^U]. Its size in the 
horizontal (x) direction is on average 120 /im. Although the beam spot size in the vertical (y) 
direction is only 5.6 ^m, we use a constraint of 50 fim in the vertex fit to account for the flight 
distance of the B° in the x-y plane. 

To suppress leptons from charm meson decays, we use only high-momentum leptons in the 
range 1.3 < pp < 2.4 GeV/c. The ir s candidates have momenta (pir s ) between 60 and 200 MeV/c. 
We reject events for which the x 2 probability of the vertex fit, TV, is less than 0.1%. We do not 
use dedicated optimization procedures for these cuts, rather we apply a selection criterion to a 
likelihood ratio, X , calculated from the signal and background distributions of pp, p 7Ts , and Vy . We 
reject events for which X is lower than 0.3 and we retain the I — tt s pair with the highest value of 
X when more than one candidate is found. 

The flavor of the other B meson (-Btag) is determined by means of lepton tagging (£ = e, /x). 
A cut on the momentum in the center-of-mass frame is applied: 1.0 < pp < 2.35 (1.1 < pp < 2.35) 
GeV/c for electrons (muons). If more than one lepton track in an event survives, the one with the 
largest momentum is used to determine the flavor of the i?tag- The track selected for tagging is 
used to compute the tag vertex position, using a procedure analogous to that described above. We 
finally require \Az\ < 3 mm and < o/\ z < 0.5 mm. 

These selection criteria accept 470,877 events in the data sample, which we will refer to as 
tagged events, and 5,291,868 partially reconstructed events that fail only the requirements of lepton 
tagging, which W6 will refer to &s untagged events, ctrid to untagged, if the B Yec is reconstructed 
as a B° and B° untagged if the B iec is reconstructed as a B°. Fig. ^ shows the distributions of 
the squared neutrino mass for tagged, B° untagged and B° untagged events. Different components 
are recognized in the total sample: signal events, due to D*~l + vp events, including the radiative 
decays D*~l + up 7/71"° and other decays with the same final signature, like B° — > D*t/X c (t/X c — > 
IX); peaking background events, that correspond to decays of charged B mesons which peak in 
the distribution, as the signal, like B + — > D** l + vp, D**° — > D*~ir + ; combinatorial events, 
corresponding to all BB decays not included in the previous categories; continuum events, coming 
from light quark decays e + e~ — > qq with q = u,d,c, s or e + e~ — ► £ + £~ with £ = e, fi, r. 

The following procedure is used to determine the fractions of the various components in data. 
The distributions are fitted separately for the different types of events, using Monte Carlo events 
for the BB components and off-resonance events for the continuum component. The probability 
density functions (PDFs) for the different components are a Gaussian plus a function B(M^) Vi\ 




(5) 
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Figure 1: Squared neutrino mass distributions for the tagged (top), B° untagged (center) and B° 
untagged (bottom) data events. The following fitted contributions are superimposed: continuum 
(solid grey line), combinatorial (dotted), B ± peaking (dash-dotted), signal (dashed), all (solid 
black). 
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for signal and peaking background, and piecewise continuous polynomials for combinatorial and 
continuum samples. For signal events, we find the distribution peaking slightly below zero mainly 
because of soft pion decays in flight and electron energy losses; as for the peaking background, the 
shift to larger positive values is due to the fact that in this case we neglect the presence of one (or 
more) additional particles. We then fit the data sample, fixing all the shape parameters for peaking, 
combinatorial and continuum background. The fraction of continuum f cont is also fixed, using 
the ratio between on-resonance and off-resonance luminosities, as well as the fraction of peaking 
background ^P eak ) that is extracted from the fit on Monte Carlo. Conversely we float the signal 
shape parameters and the fraction of signal / Slg , determining also the fraction of combinatorial 
background / comb = 1 - / sig - / peak - / cont . 

The fractions determined with this procedure are summarized in Table ^ for the different data 
samples. The M 2 region —4.0 < M 2 < 2.0 GeV 2 /c 4 is defined as signal region and will be used in 
the nominal fits from now on. The remaining part (-10.0 < M 2 < -4.0 GeV 2 /c 4 ) will be referred 
to as side-band region and will be used in some fits to At to determine combinatorial background 
parameters. Tabled also shows fractions extrapolated to the signal region only. 

Table 1: Fractions of signal, continuum, peaking and combinatorial background, extracted from 
My fits for the tagged, B° untagged and B° untagged samples. The first uncertainties refer to 
the statistical uncertainty from the fit, the second are systematic errors from the shape parameter 
fixing. 



Sample 


Parameter 


Fraction in the 


Fraction extrapolated 






whole M 2 region 


to the signal region 


Tagged 


/Slg 


(32.72 ± 0.15 ± 0.58)% 


(46.59 ± 0.21 ± 0.82)% 




jpe&k 


4.03% 


5.54% 




y^cont 


4.12% 


3.14% 




ycomb 


(59.13 ± 0.15 ± 0.58)% 


(44.73 ± 0.21 ± 0.82)% 


Untagged + 


/S1 g 


(23.51 ± 0.03 ± 0.27)% 


(35.72 ± 0.05 ± 0.41)% 




jpe&k 


3.05% 


4.47% 




j cont 


6.47% 


5.38% 




ycomb 


(66.97 ± 0.03 ± 0.27)% 


(54.43 ± 0.05 ± 0.41)% 


Untagged — 


/S1 g 


(24.23 ± 0.04 ± 0.20)% 


(36.33 ± 0.06 ± 0.30)% 




ypeak 


3.02% 


4.42% 




y^cont 


6.14% 


5.20% 




ycomb 


(66.61 ± 0.04 ± 0.20)% 


(54.05 ± 0.06 ± 0.30)% 



4 TIME-DEPENDENT ANALYSIS 

Since the B° mesons can undergo oscillation, we will refer to mixed events when both decay in the 
same flavor state (B°B° or B°B°), and define the mixing state s m = — 1 for such events. Conversely, 
we will refer to unmixed events (s m = 1) if the flavors of the two B mesons are opposite. Similarly 
we define St = 1 if the -B tag is a B°, st = — 1 if the 5 tag is a B°. 

The calculation of the decay rate for the four tag-mixing states is straightforward, since we are 
dealing with semileptonic B° decays both on the reconstructed and tag side and these are pure 
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flavor eigenstates: 



jr(At) 



2t b o 

I I At 



-2s* 



[1 - cos(Am d At)] 



s m = 1 : F(At) = e r s° [i + C os(Am d At)] (6) 

2t b o 

where Amj is the B° oscillation frequency, t b o is the average lifetime of the physical states and we 
set Ar, the lifetime difference between the physical states, equal to 0. 

For the At distribution of signal events we use Eq. ©, modified to account for several experi- 
mental effects. 

The -B-tagging algorithm introduces a modification to the theoretical PDFs as there is a finite 
probability per event (called the mistag rate) of falsely tagging with a wrong-sign lepton candidate. 
Considering a B° mistag rate (i.e. the probability for a true B° being tagged B°) and a B° 
mistag rate (i.e. the probability for a true B° being tagged B°) , we define w as the average value 
of these probabilities and Aw as their difference. Also, the signal PDFs are modified to reflect the 
fact that reconstruction efficiency for ^ + vr~ pairs can be different from the reconstruction efficiency 
for £~ir^ pairs. We therefore define an average reconstruction efficiency s rec and a reconstruction 
asymmetry A rec . In the same way, the tagging efficiency for positive leptons can be different from 
the tagging efficiency for negative leptons, so we define an average tagging efficiency Et ag and a 
tagging asymmetry Ai ag . We normalize our distributions to an average reconstruction efficiency 

£-rec — 1- 

To account for the finite resolution of vertex determinations, the function of the true lifetime 
difference Att rU e must be convolved with a decay time difference resolution function. We adopt 
a 3-Gaussian description where the two main contributions (referred to as "narrow" and "wide") 
depend on the per-event error of the vertex separation, a At, while the third ("outlier" Gaussian) 
is independent of o~At'- 

1 

T^ipn, Snibwi Swi So, fwi fo) = fo 7== 6 2s ° 4" 

y/ZTTSo 



1 (t-hm<T Ai ) 2 -i (t-ft n tT At ) 2 

+f w ^= e + (1 - f w - f )-= e (7) 

V 2irs w a At V 2irs n a At 

where t = At — At true , bi and Si {i = n,w) are respectively the biases and the error scale factors 
for the narrow and wide components, while s Q is the standard deviation of the outlier Gaussian 
(for this component we assume b a = 0) ; f w and f Q are respectively the fractions of wide and outlier 
Gaussians. 

Accounting for these effects, we define a PDF for signal events and direct leptons, meaning 
tagging leptons that do not originate from secondary charmed meson decay: 

^(At,a At \s m ,s t ) = F meas {At,o At \s m ,s t )®K% (8) 

where, using the definition k = \q/p\ — 1: 

s m = -l: ^ mcas (At) = i^_(l + s t A rec )e'^ {[1 - s t (Aw + 2k(l - w) - A tag (l - 2w))\ 

4r B o 

-[l-2w + st(A ta g - 2k(l - w))] cos(Am d At)} 

(9) 
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s m = l: F mc ™(At) = i*L(l - s t A rec )e {[1 - s t (Aw - 2kw - A tag (l - 2w))) 

At b o 

+ [1 - 2w + s t (A tag + 2kw)} cos(Am d At)} (10) 

For signal events in which the tagging lepton comes from a secondary charm decay (cascade 
leptons), we use different PDFs according to the two possibilities that the charm meson is a product 
of the -Btag or is the unreconstructed D° on the decay side (decay-side tagged events). In the 
former case we just take the signal PDFs with the reversed tag information: ^^(At, a A t\s m , St) = 
J-V* (At, a A t\ — s m , —St) and we allow for different mistag and resolution parameters. In the latter 
case, we use a simple exponential term: 

X _ I A *l 

^bta g (Ai, a A t\s m , s t ) = - e r o e (l + s m L> Dtag )(l + s t A' t )(1 - s t s m A rec ) ® TZ Dtag (11) 

where T£> e is an effective D° lifetime and A' tag is a distinct tagging asymmetry for this sample (see 
Sec. EJ). We characterize from now on the mistag in a given sample j in terms of the dilution 
parameter Dj = 1 — 2wj. 

For charged B peaking background events, we similarly use a pure lifetime PDF for the direct 
component: 

1 _JM. 

^ ir (At, <T At \s m , s t ) = e t b+ (1 + s m D + )(l + s t A tag )(l - s t s m A rec ) <g> (12) 

2t b + 

while we use Eq. ((TT|) for the decay-side tagged events. 

For combinatorial background events, we recognize different components that can be parameter- 
ized with the PDFs defined above, but with sets of effective parameters that need to be determined 
from the sidebands. For the neutral part we take into account three contributions: first, the 
events in which both leptons in the events originate directly from the two B mesons, even if the 
reconstructed side is not a genuine D*~£ + vi event. For example, these may be other semilep- 
tonic B decays in which the reconstructed soft pion candidate is actually a random track in the 
event. For these events we use Eqs. © and (|10|) . where some shape parameters are allowed to 
differ from the corresponding ones for signal, since the vertex position is mis-determined using a 
wrong pion track. We also allow Area to be distinct from A^ c for the same reason, but we use the 
same CP asymmetry parameterization of the signal, since this kind of events carries the correct 
tag-mixing information. Second, the same procedure applies to cascade combinatorial events. The 
third contribution comes from events in which the tag track is coming directly from a B, while the 
lepton on the decay side comes from a secondary charm meson. For this kind of events we assume: 
•^cas g 2(^' a At\sm, St) = ( At, dAt I — Sm, St) and we allow for different mistag and resolution 
parameters with respect to the direct lepton case. For decay-side tagged and charged B events 
the descriptions of Eqs. (|11|) and Q12j) apply also for combinatorial events, given the changes of 
parameterization for resolution and asymmetries. 

For continuum events we use a single lifetime distribution: 

^ cont (At, a At \s m , s t ) = e~w (l + s m D cont )(l + s t A%f )(1 - s t s m A c r °?) ® K^ r s (13) 

•^Tcont 

Its parameters are all determined from a fit to off-resonance events. 
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In summary, the At total PDF will be a sum of all the terms so far introduced: 



.F(Ai,cr At |s TO ,s t ) = 

/ SiS {(l " «&,)[(! " fl&)*& + ^cal^af] + fl&^W 
+/ Peak [(l-5S?a g )^ r +^fa g ^"Dta g ] 

+ / C ° mb {(l " ^ S " $g*) [(1 " ~ 9^)^ 

i _bkgrbkg i bkg ^-bkg -i bkg ^-bkg . bkg ~- \ 
"Ti/cas ^cas ~r i/ C as,2- r cas,2J 9+ • r + i/Dtag^Dtag/ 

_|_ycont^-cont (14) 

where <?j are generic fractions of the tagging category j in sample i and we assume the fractions of 
decay-side tagged events, determined from Monte Carlo simulation, to be equal in the signal and 
peaking background components. 

A crucial point of this analysis is to determine particle detection asymmetries, that could in 
principle be degenerate with the CP asymmetry. A method commonly used for determining detec- 
tion asymmetries in time-dependent analyses is to use tagged and untagged events reconstructed 
as B° or B°, to have four equations from which etag, e rec , A rec and At ag are determined uniquely. 
This procedure, fully explained in Appendix of [7], is however appropriate when CP violation in 
mixing is assumed to be zero. Here we must employ a different approach. We parameterize the 
total number of expected events for each tagged and untagged category in a likelihood function, 
which contains a term coming from the shapes of the PDFs for tagged events, an extended term for 
the number of tagged events and an extended term for the number of untagged events. We show 
the procedure in full detail for signal events. The expected number of tagged events per category 
is obtained in a straightforward way by integrating Eqs. ® and (jl()j) . while for untagged events: 

N exp (s u ) = N tot - S " ^ [l - e tag + S u EtagA tag X d - S u k(l - X d )(l - £ ta g)\ (15) 

where Ntot is the total number of events (tagged and untagged), xj, = 1/[1 + (r B oAmrf) 2 ] and 
s u = +1(— 1) for B° (B>°) untagged events. 

The likelihood function to be maximized is of the form: 

\nL tot = \nL- £ In L ext (s t , s m ) - £ In L ext (s u ), (16) 

where \nL ext {i, (j)) = ln(N^p\) — N^' InNlxp' 1 + Nl^p are the extended terms for tagged and 
untagged events while InL is the log-likelihood constructed starting from Eq. (|14jl. This method 
reduces the correlation between CP mixing asymmetry and particle detection asymmetry, that 
nevertheless remains significant, as quoted in Sec. El 



5 FITS OF SIMULATED SAMPLES 

Several tests were performed on simulated events to verify that the fit determines the free parame- 
ters k, A rec and At ag correctly. First, fits to the various components of the Monte Carlo sample and 
to the total BB sample were performed to check for the presence of biases in the analysis technique. 
This Monte Carlo sample was generated with a value of k equal to 0. We summarize in Table |^1 
the result of these fits for side-band region (first column) and signal region (second column), r^o, 
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t b + and AiTij are fixed to the generated values in this procedure, while the normalizations of the 
various contributions are fixed to the values extracted from the distribution. 

The most notable parameters in Tableware the asymmetries k, A rec , Atag, Arec and A[ ag . We 
find all to be compatible with 0, except Arec and A' tag . The source of Arec 7^ is random tracks 
reconstructed as soft pions in the combinatorial background. We do not use a control sample for 
estimating Arec but we use the value found in the side-band fit; its possible variation from side-band 
to signal region is taken into account during systematic evaluation. From Monte Carlo studies, we 
find that the source of A' tag 7^ is charged kaons from the unreconstructed D° faking muons. We 
compare the asymmetry value found in the Monte Carlo sample (second column of Table |2| with 
the value calculated in a charged kaon control sample [D* — > D°tt, D° — > Kir decays) from the 
integrated K + and K~ efficiencies over the tag lepton momentum spectrum. This is found to be 
compatible with the fitted value (A' tag cs = (18.0 ± 0.5)%). We then conclude that the asymmetry 
can be considered as entirely due to this source and we fix the value of A' tag to the one found in 
the control sample. A similar procedure is then exploited in the data fit. 

Other checks are performed using large numbers of toy experiments, each generated with statis- 
tics corresponding to the size of the data sample. A set of 120 toy experiments is used to test the 
validity of the uncertainty calculated by the minimization algorithm and to check that the final 
value of the negative-log-likelihood is compatible with the value found in the data. Other Monte 
Carlo samples are generated with values of k or A rec or At ag equal to ±0.01 to check the feasibility 
of disentangling the asymmetries from one another. The fits show that this disentanglement is 
possible. Also, a CP asymmetry different from is generated in the generic B°B° sample to verify 
the sensitivity of the fit to values of k 7^ 0. We find k = (8.9 ± 2.9) • 10 -3 for a generated value of 
10 • 10~ 3 and k = (-10.1 ± 2.9) • 10~ 3 for a generated value of -10 • 10~ 3 . 

6 FITS OF THE DATA SAMPLE 

The same procedure used in Monte Carlo events is applied to the on-resonance data. We float 
the same set of parameters and we fix the others in the following way. The tagging efficiency £t ag 
is determined from the number of tagged events divided by the total number of events: we find 
£tag = (9.06±0.02)%. For continuum events we use parameters taken from the off-resonance fit. For 
combinatorial background, some parameters are floated, some are fixed from a previous side-band 
fit (see Table |2J). For signal and peaking background, we take dilutions from Monte Carlo and we 
float most resolution parameters. For A' tag , we approximate its value using the estimate from the 
control sample explained above, assuming the fraction of charged kaon mistagging is the same as 
in Monte Carlo (we tested this hypothesis, as described in Sec. [7J. We obtain an asymmetry of 
A' tag = (7.6 ±0.8)% from the control sample. t b o, Amj and t b + are fixed to their PDG values [Tl] : 
the former two are then floated to estimate the fit stability in data when these parameters are free. 

The fit results are shown in Table (third and fourth column); for the three signal asymmetry 
parameters we find the following correlation coefficients: p(Af ag , A re c) = 43.4%, p(k, A re c ) = 50.2%, 
p(k,A ta g) = 88.5%. Figs. 121 and El show the fitted PDFs for the subsamples equivalent to the four 
tag-mixing states; Pearson's test gives a x 2 °f 682 for 385 degrees of freedom. Fig. 0] shows the 
PDF asymmetry derived from the fit in Fig. [2] between mixed events with St = 1 and St = —1. 
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Table 2: Top: results of the fits to the BB Monte Carlo and for: side-band region (first column), 
signal region (second column) . Results of the fits to the on-resonance sample for side-band region 
(third column), signal region (fourth column). We also report asymmetrical errors for some fun- 
damental parameters. A horizontal line is used to separate signal from background parameters. 
Bottom: Continuum parameters determined from the total off-resonance sample: they are fixed in 
the nominal on-peak data fit. 



Parameter 


Side-band 


Signal region 


Side-band 


Signal region 






fit (MC) 


fit (MC) 


fit (data) 


fit (data) 




k 


CO 8 ± 3 OVlO -3 


(-1 6 ± 1 7V10~ 3 


(-2 1+4 9V10~ 3 


(6 ± 3 4V10~ 3 [+3 4 -3 
(1.0 ± 1.5)-10~ 3 [+1.5, -1 
(3.1 ± 2.1)-10~ 3 [+2.1, -2 


51 


A 




(0.02 ± 0.06)-10~ 3 




51 


Atag 


(-0.2 ± 2.4)-10~ 3 


(0.3 ± 1.1)-10~ 3 


(-0.2 ± 2.9)-10~ 3 


11 


b dir 




-0.049 ± 0.002 




-0.044 ± 0.008 








0.859 ± 0.007 




1.15 ± 0.02 




fw 




(9.6 ± 0.3)% 




(9.2 ± 0.7)% 




Sw 




1.90 ± 0.03 




1.90 (fixed) 




So 


_ 


(7.2 ± 0.5) ps 


_ 


7.2 ps (fixed) 




So 


- 


(2.7 ± 0.3)-10~ 3 


- 


2.7-10^ 3 (fixed) 




Wdir 


- 


(1.32 ± 0.04)% 
(-3.2 ± 7.8)-10~ 4 


- 


1.32% (fixed) 
-3.2-10" 4 (fixed) 




AWdir 


- 


- 




TD C 


(0.226 ± 0.009) ps 


(0.267 ± 0.008) ps 


(0.342 ± 0.015) ps 


(0.321 ± 0.010) ps 




^cas 


- 


(20.8 ± 2.4)% 


- 


20.8% (fixed) 






- 


(1.0 ± 0.4)% 


- 


1.0% (fixed) 




teas 
u n 


- 


-0.33 ± 0.04 


- 


-0.28 ± 0.09 




icas 


- 


-2.4 ± 0.3 


- 


-2.4 (fixed) 




A ta g 


- 


(18.9 ± 1.0)% 


- 


7.6% (fixed, see text) 




-^Dtag 


- 


(83.7 ± 2.0)% 


- 


83.7 (fixed) 




yeas 


- 


(6.67 ± 0.20)% 


- 


(5.0 ± 0.5)% 




#Dtag 


- 


5.62% (fixed) 


- 


5.62% (fixed) 




D+ 


(80.6 ± 0.8)% 


(90.1 ± 2.5)% 


(83.2 ± 1.8)% 


(83.4 ± 4.2)% 




bk g 

T 


(1.469 ± 0.017) ps 


(1.28 ± 0.04) ps 


(1.51 ± 0.03) ps 


(1.53 ± 0.02) ps 




A bkg 

Am d 


(0.479 ± 0.007) ps" 1 


(0.37 ± 0.04) ps" 1 


(0.506 ± 0.010) ps" 1 


(0.586 ± 0.010) ps" 1 




4 bkg 


(3.9 ± 2.3)-10~ 3 


3.9- 10~ 3 (fixed) 


(9.7 ± 2.9)-10~ 3 


9.7-10" 3 (fixed) 




»S? 


(1.19 ± 0.07)% 


1.19% (fixed) 


1.19% (fixed) 


1.19% (fixed) 






(1.4 ± 0.9)-10~ 3 


1.4-10~ 3 (fixed) 


1.4-10" 3 (fixed) 


1.4-10~ 3 (fixed) 




7 dir 
°bkg 


0.071 ± 0.036 


-0.034 ± 0.030 


0.013 ± 0.006 


-0.056 ± 0.008 




dir 
s n,bkg 


0.837 ± 0.006 


0.837 (fixed) 


0.850 ± 0.011 


0.850 (fixed) 




fbkg 
J w 


(7.6 ± 0.5)% 


(6.1 ± 0.5)% 


(12.7 ± 0.8)% 


(7.4 ± 1.5)% 




r bkg 


0.205 ± 0.009 


0.205 (fixed) 


0.420 ± 0.012 


0.420 (fixed) 




w bkg 


(18.7 ± 0.6)% 


18.7% (fixed) 


18.7% (fixed) 


18.7% (fixed) 




A™ c b a k s g 


0.017 ± 0.006 


0.017 (fixed) 


0.017 (fixed) 


0.017 (fixed) 




icas 
°n,bkg 


-0.02 ± 0.06 


-0.02 (fixed) 


-0.13 ± 0.08 


-0.13 (fixed) 




i cas,2 
°n bkg 


0.02 ± 0.08 


0.02 (fixed) 


0.15 ± 0.11 


0.15 (fixed) 




bkg 
i/Dtag 


(13.0 ± 0.2)% 


(9.2 ± 0.8)% 


(8.5 ± 0.5)% 


(6.1 ± 0.8)% 






(40.7 ± 0.3)% 


(37.1 ± 0.2)% 


(33.7 ± 0.5)% 


(31.9 ± 0.3)% 





Parameter Value from 
off-peak fit 
(0.62 ± 0.03) ps 
A c r °f (6.4 ± 2.1)% 

A c ta f (3.9 ± 2.9)% 

Dcont (47.8 ± 1.6)% 
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Figure 2: Fit to the At distributions for the on-resonance data samples in linear scale. Letters 
indicate the different samples: (a) Unmixed positive (s m = l,st = 1), (b) Unmixed negative 
(s m = I, St = —1), (c) Mixed positive (s m = — l,St = 1), (d) Mixed negative (s m = —l,st = — 1). 
The following fitted contributions are shown in the fit: continuum (dark grey), combinatorial (light 
grey), peaking (black), signal (white). 




A t (PS) A t (ps) 



Figure 3: Same as Fig. [21 in logarithmic scale. 
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Fi gure 4: Raw asymmetry between mixed positive (s m — — ljSt — 1) and mixed negative (s m — 
—l,st = —1) events with the PDF asymmetry derived from the fit in Fig. 21 superimposed. The 
asymmetry at low values of \At\ corresponds to the kaon-mistag contribution in decay-side tagged 
events (see text). 
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7 SYSTEMATIC STUDIES 

We consider the following sources of systematic uncertainties. 

For the reconstruction asymmetry for combinatorial background events, we repeat the fit in 
the first column of Table 21 for combinatorial Monte Carlo events in the signal region, finding 
A re c = (6.5±2.1) • 10 -3 . The difference between this result and the one in the sideband (= 2.6- 1(T 3 ) 
is applied to the reconstruction asymmetry in data and the fit is re-done to estimate the impact 
on k. 

We vary the continuum reconstruction and tagging asymmetries by ±lcr, where the standard 
deviations are taken from the fit to the off-peak events. The variations in k are averaged and added 
in quadrature to give the systematic uncertainty. 

We perform a fit in which the asymmetry A' tag is free rather than fixed from a control sample: 
we obtain A' tag = (4.0 ± 3.8)% which is compatible with the assumption that the fraction of 
kaon mistagging is the same as in Monte Carlo, but the uncertainties are large because of the high 
correlation between At ag and A' tag . The resulting variation in k is taken as a systematic uncertainty. 

We vary the sample fractions by ±1<7, where the statistical errors are taken from the data 
fits, both to tagged and untagged events separately. Moreover, for B° and B° untagged events, we 
move the two central values in opposite directions in order to consider the maximum asymmetry 
that could be generated by the fit uncertainty. Similarly, we vary the fractions of their correlated 
systematic uncertainties simultaneously for tagged and untagged events. Continuum fractions are 
varied by the uncertainty in the on/off-resonance luminosity ratio, which we take as 1.3% from |lUj . 

For the analysis technique, we determine the mean of k^ t — k gen in the Monte Carlo tests we 
performed and we correct the result on data accordingly; the full size of the correction (5.5 • 10~ 4 ) 
is taken systematic uncertainty. 

We let t b o and Am^ float in the data fit, obtaining t b o = (1.551 ± 0.012) ps and Am,/ = 
(0.475 ± 0.019) ps" 1 which are in good agreement with the world averages [Hj, considering a 
known bias on Am^ caused by the partial reconstruction technique |15j . 
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Most mistag parameters for the various samples are fixed in the data fit. When omitting 
decay-side tagged events, we can assume mistagging is almost entirely due to charged pions. Using 
inclusive pion control samples, we notice that for both negative and positive pions, the relative 
difference between data and Monte Carlo is smaller than 10%. We conservatively vary the mistags 
of 10% and quote the variations of A: as a systematic uncertainty. 

Also some parameters of the resolution are fixed, namely the fraction and width of the "outlier" 
Gaussian and the width of the "wide" Gaussian are fixed. For the former, we repeat the fit setting 
jo = (2-Gaussian resolution only), the latter is varied by 30% conservatively. 

All other parameters that are fixed from the sidebands or from the off-resonance fit are varied 
by their statistical errors to estimate corresponding systematic uncertainties. 

All the variations in k are summed in quadrature. The effect on the parameter k for all these 
sources, as well as the total uncertainty is summarized in Table 03 



Table 3: Sources of systematic uncertainty for the \q/p\ measurement. 



Source Syst. error (-10 -3 ) 

Reconstruction asymmetry for combinatorial background 1.1 

Asymmetries for continuum 1.0 

Tagging asymmetry for decay-side tagged events 0.2 

Ml fractions 1.0 

Likelihood fit bias 0.6 

Physical parameter fixing 0.5 

Mistag parameters 0.0 

Resolution function 0.3 

Sideband-fixed parameters 0.5 

Total 2l) 



8 RESULTS 

From the fit performed on the on-resonance data sample we measure: 

\q/p\ - 1 = (6.5 ± 3.4(stat.) ± 2.0(syst.)) • 10" 3 

which relates to CP violation in mixing. 

To facilitate comparison with other measurements and formalisms, we also express our result 
in two other commonly used notations. The parameter A$l defined in Eq. Q is found to be: 

1 _ | /r>| 4 

A SL = 1 - ^ 4 = (-13.0 ± 6.8(stat.) ± 4.0(syst.)) • lO" 3 
or, using the eb parameter and the relation q/p = (1 — £b)/(1 + £b)'- 

Ree B ^^SL = ( _ 3>2 ± L7 ( stat .) ± i. (syst.)) • 10~ 3 



1 + \£B 
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9 SUMMARY 



We have presented a determination of the parameter \q/p\, that is a measurement of CP violation 
m B°B° mixing, using a fit to At, the time difference between the two B decays. One of the B 
mesons is partially reconstructed in the semileptonic channel D*~£ + V£, i.e. only the lepton and the 
soft pion from D*~ — > D°tt~ decay are reconstructed, while the flavor of the other B is determined 
by means of lepton tagging. We use a luminosity of 200.8 fb _1 collected by the BABAR detector in 
the period 1999-2004, and obtain the preliminary result: 

\q/p\ - 1 = (6.5 ± 3.4(stat.) ± 2.0(syst.)) • 10~ 3 . 

This result is compatible with the current world average, and the magnitude of its error is 
comparable to those of the most recent other measurements [Sj- The corresponding central 
confidence interval at 95% C.L. for \q/p\ is [1.0012 - 1.0142]. 

10 ACKNOWLEDGMENTS 

We are grateful for the extraordinary contributions of our PEP-II colleagues in achieving the 
excellent luminosity and machine conditions that have made this work possible. The success of 
this project also relies critically on the expertise and dedication of the computing organizations 
that support BABAR. The collaborating institutions wish to thank SLAC for its support and the 
kind hospitality extended to them. This work is supported by the US Department of Energy and 
National Science Foundation, the Natural Sciences and Engineering Research Council (Canada), 
Institute of High Energy Physics (China), the Commissariat a PEnergie Atomique and Institut 
National de Physique Nucleaire et de Physique des Particules (France), the Bundesministerium fur 
Bildung und Forschung and Deutsche Forschungsgemeinschaft (Germany), the Istituto Nazionale di 
Fisica Nucleare (Italy), the Foundation for Fundamental Research on Matter (The Netherlands), the 
Research Council of Norway, the Ministry of Science and Technology of the Russian Federation, and 
the Particle Physics and Astronomy Research Council (United Kingdom) . Individuals have received 
support from the Marie-Curie IEF program (European Union) and the A. P. Sloan Foundation. 

References 

[1] M. Ciuchini et al, JHEP 0308, 031 (2002). 

[2] M. Beneke et al, Phys. Lett. B 576, 173 (2003). 

[3] S. Laplace et al, Phys. Rev. D65, 094040 (2002). 

[4] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

[5] K. Abe et al, The BELLE Collaboration, Phys. Rev. D73, 112002 (2006). 

[6] B. Aubert et al, The BABAR Collaboration, Phys. Rev. Lett. 96, 251802 (2006). 

[7] B. Aubert et al, The BABAR Collaboration, Phys. Rev. D70, 012007 (2004). 

[8] |http://www-d0.fnal.gov/Run2Physics/WWW/results/prelim/B/B29/B29.pdf| 

21 



[9] M. Bona et al, The VTfit Collaboration, hep-ph/0605213 

[10] The BABAR Collaboration, B. Aubert et al, Nucl. Instrum. Methods A479, 1-116 (2002). 

[11] G. C. Fox and S. Wolfram, Phys. Rev. Lett. 41, 1581 (1978). 

[12] Unless specified, all four-momenta are measured in the T(4S) rest frame. 

[13] With m = M*, we define B{m) = [{m - m ) ai + c(m - m ) a2 ]e- b ( m - m °\ 

[14] Particle Data Group, S. Eidelman et al, Phys. Lett. B592, 1 (2004). 

[15] B. Aubert et al, The BABAR Collaboration, Phys. Rev. D73, 012004 (2006). 



22 



